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Abstract 

We show that graphs, networks and other related discrete model sys- 
tems carry a natural supersymmetric structure, which, apart from its con- 
ceptual importance as to possible physical applications, allows to derive a 
series of spectral properties for a class of graph operators which typically 
encode relevant graph characteristics. 



1 Introduction 



In the following we show that graphs and networks support in a very natural 
way a supersymmetric structure. Our starting point is a class of geometrically 
relevant graph- operators acting on the direct sum, 7i, of the vertex-Hilbert space, 
Ho and the edge-Hilbert space, TCi- As to these technical prerequisites cf. the 
analysis being made in pQ and further references given there. 

Of particular importance is a graph-Dirac operator, D, introduced there, 
which maps Tio into Hi and vice versa. Furthermore we have a natural Laplace 
operator on H which together with D and another supercharge forms a closed 
(N = 2)-superalgebra. This abstract susy-structure, if concretely represented 
by our graph-operators, allows, among other things, to derive a series of useful 
spectral properties of these operators. 

The natural existence of this susy structure on graphs (and related models) 
may also be of some relevance in a wider context. We have been promoting a 
discrete network approach to quantum space-time physics in recent years (see 
for example jT^j or ^1 and further references given there) which is assumed 
to underly our more macroscopic continuum physics on the Planck-scale. Our 
present analysis shows that these discrete model systems, perhaps contrary to 
naive wisdom, are in fact quite rich as to their structural properties. 

We conclude this introduction with the remark that such technical properties 
of graphs have also been mentioned elsewhere in a however different context and 
using a different formalism (cf. [U|). So we should add the remark that there 
may exist papers in other fields of research being related to our work but which 
we presently are not aware of. 

2 Some Concepts from Algebraic Graph Theory 

To set the stage, we briefly introduce some concepts employed in algebraic graph 
theory but in a form particularly adapted to our own framework (see also pQ; 
as to the wider context cf. [2] to jH] and further references given there). The 
graph, G, consists of a countable set of vertices (or nodes), V, with Xi the labelled 
vertices, and an edge set, E, the directed edges, pointing from Xj to Xj, denoted 
by dij. We assume (for convenience) the graph to be simple, that is, an edge, dij, 
can be associated with the ordered pair of vertices, (xi,Xj) with i ^ j. 

So-called unoriented graphs with the edges consisting of unordered pairs, 
{xi,Xj}, can be subsumed in the above framework in two different ways, each 
of which having a certain advantage of its own. On the one hand, we can 
give the graph an arbitrary orientation, that is we associate to each undirected 
edge,{xj, Xj}, one of the two possible choices, (xi, Xj) or (xj, Xi). It turns out that 
most of the concepts and calculations do not depend on the particular choice (see 
below). On the other hand, we can associate an unoriented but orientable graph 
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with a directed graph so that each edge occurs twice, that is if dij belongs to E, 
dji also belongs to E. 

For reasons of simplicity we assume our graph to be locally finite, that is, 
each vertex is only incident with a finite number of edges. A slightly stronger 
assumption is it to assume this vertex degree to be globally finite over the vertex 
set. For a directed graph we then have ingoing edges and outgoing edges relative 
to a given node with the respective vertex degrees, v * n and v° ut and the total 
degree v t = v\ n + v° ut . 

For such a graph we can introduce two Hilbert spaces, a vertex Hilbert space, 
Ho, and an edge Hilbert space, Hi, with orthonormal bases the set of vertices, 
Xi, and the set of directed edges, d^. This means, we introduce a formal scalar 
product on Ho, Hi respectively with 

{Xi,Xj) §ij , {dij, dim) $il$jm iX) 

and with vectors being the formal sums 

00 oo 

f = ^2fiXi , 9=^2 9n d H with /*' 9ij e C (2) 

1 i,j=l 

with Yl \fi\ 2 < 00 and Yl \9ik\ 2 < oo. 

Remark: We treat the vertices and edges as abstract basis elements (in a way 
similar to the group algebra of a group). One can of course consider the abstract 
vectors equally well as discrete functions over the vertex- or edge-set, respectively 
and the basis vectors as elementary indicator functions. 

If we deal with an undirected but orientable graph, employing the second 
variant introduced above, we found it convenient (cf. ilj) to introduce the super- 
position 

bij . d{j dji bji (3) 
and relate it to an undirected but orientable edge. 

3 Operators on Graphs 

In [T] we introduced two operators, interpolating between Ho and H\. We define 
them on the basis vectors: 

d{xi) := <ki - d ik' (4) 

k k' 

with the first sum running over the ingoing edges relative to x%, the second sum 
running over the outgoing edges. In the case of a symmetric (or undirected graph; 
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version two) we have 



d{xi) := ^ ( d ki ~ dik) = ^2 bki ( 5 ) 



A- 



This operator is closely related to a sort of non-commutative discrete differential 
calculus on graphs as we have 

df = - m k (6) 

i,k 

A simple calculation shows that the adjoint, d*, acts on the basis vectors of Hi 
as follows: 

d*(d ik ) = x k - Xi (7) 

In algebraic graph theory (finite graphs) the so-called incidence matrix, B, 
is introduced, having the entry 1 if vertex Xi is the positive end of a certain 
(ingoing) edge, and having a —1 if it is the negative end (outgoing edge) (see for 
example 3J). This matrix corresponds to our operator d*. 

Another important operator is the adjacency matrix, A, being a map from H.Q 
to 7io and having (in ordinary graph theory of (un)oriented graphs) a +1 at entry 
if Xi and Xj are connected by an edge. This matrix is a symmetric operator, 
In our more general context (which includes however the ordinary 
situation as a special case) of directed graphs one can introduce the in-adjacency 
matrix, A in , and the out- adjacency matrix, A out , with A = A in + A out . In our 
(operator)-notation they are given by 

Axi = ^ ] t k j x k , A Xi = ^ ^ Xk , A Xi = ^ ] x^ (8) 

k— *i i—*k' 

with ~ designating the unordered pair {xi,Xk}, k — > i the ordered pair (k,i) 
and e^i is either one or two depending on the two possible cases of one directed 
edge between node Xi and node Xk or two directed edges, pointing in opposite 
directions. 

These operators can be built up from more elementary operators (cf. pQ). 

d\Xi = ^2 dki , d 2 Xi = 22 dik' (9) 



so that 



d\ d ik = Xk , d* 2 dik = Xi (10) 
d = d\ — di , d* = d\ — d* 2 (11) 
d\d x Xi = vf 1 ■ Xi , d* 2 d 2 Xi = v° ut ■ x { (12) 
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d\d 2 Xi=y] x k > , d* 2 d 1 x i = J ^2x k (13) 



i— >fc' k— *i 

An „,out 



where v™ , v° is the in-, out degree of vertex x^ respectively. We hence have 
Lemma 3.1 The in-, out-vertex degree matrices read 

V in = d\d l , V out = d* 2 d 2 (14) 
The in-, out-adjacency matrices read 

A in = d* 2 d 1 , A out = d\d 2 (15) 
A = A m + A out is symmetric. 

Proposition 3.2 The so-called graph Laplacian is the following positive operator 

-A := d*d = (V in + V out ) - {A in + A out ) — V — A (16) 

The reason to call this operator a Laplacian stems from the observation that 
it acts like a second order partial difference operator on functions of H . 



i \ k^>i i^>k / 



(17) 



and after a simple relabelling of indices 

- a / = - £ ( j2 h + E /* - - «r7< ) * 



\k—*i i— >k 



- E ( E(/* - /*) + Eta - /<) ) ^ = - E ( E ^(/* - /<) ) 



x» (18) 



which reduces to the ordinary expression in the undirected case. 

Forming now the direct some TC := TCo&)'Hi, we can introduce yet another im- 
portant graph operator which closely entangles geometric and functional analytic 
properties of graphs (and similar structures); see pQ. 

Definition 3.3 We define the graph Dirac operator as follows 

d* \ TT {Ho 



D-.n^H withD := y d Q j , H=y£) ( 1 9) 

Observation 3.4 

D2 = DD =( d l d £r\ 



with d*d = —A. 
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The action of dd* on a basis vector dik reads 

dd* d ik = d (x k - x{) = ^ d k>k - ^2 dkk " ~ ^2 dli + diV ( 21 ) 

k' k" I I' 

which, after some relabelling and introduction of the Kronecker delta function 
can be written as 

dd d^ ^ {dmj&jk dj m 8jk d m j8ij -\- dj m 8ij^ (22) 

m,j 

For a function g = 9ikdik we hence get 

km (23) 



dd * 9 = ^2 I ^2 9im ~ 9mi ~ 9il + 9li ) dh 

l,m \ i / 



In the simple case of the one-dimensional directed lattice, Zi, with directed 
edges only pointing from X{ — > x i+ i, we get a transparent expression: 

dd* 9i,i+ldi,i+l \ = — ~^2 (di+l,i+2 + 9%-l,i ~ ^9i,i+l) d>i,i+l (24) 

i.e., it represents also a discrete second derivative operator on the level of edges. 
For the directed lattice, Z n , it can be related to what is called the vector Laplacian 
in the continuum (cf. [7j. p. 575). 

There is a pendant in the calculus of differential forms on general Riemannian 
manifolds where, with the help of the Hodge-star operation, we can construct a 
dual, 8, to the ordinary exterior derivative. The generalized Laplacian then reads 

-A = 8d + d8 (25) 

with 8 (modulo certain combinatorial prefactors) corresponding to our d* (see for 
example t 8j or P). 



4 Susy on Graphs 

We introduce the following simple model of a supersymmetric algebra, consisting 
of two hermitean charges, Q\ , Q 2 , and a Hamiltonian 

H s = Ql = Ql hence [H s , Q 1}2 \ = (26) 
{Qi,Qz} = hence {Q h Qj} = 2H S ■ 8^ (27) 

Defining Q + , Q_ as 

Q + :=2-\Q 1 +iQ 2 ) , Q_ := 2~ 1 {Q 1 - iQ 2 ) (28) 
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implying 

Qi = Q+ + Q- , Q 2 = -i(Q+-Q-) (29) 

we get 

Ql = Q 2 _ = , H S = {Q + ,Q„} , [H 8 ,Q±] = (30) 

That is, these three generators create a closed (N=2) susy-algebra. 

We now make the following correspondence with our graph operators: 

H s = EC = DD=^ d (31) 

yielding 



In other words, the charge Qi is our original Dirac operator. We can now check 
all the above abstract (anti) commutation relations and find that they are fulfilled 
by our representation in form of graph operators. We see that Q + is essentially a 
map from the subspace TLq to the subspace Ti\ and Q- from Hi to TCo while Q + 
vanishes on TCi, Q~ on 7Yo- Therefore we can tentatively associate TLq with the 
bosonic and TCi with the fermionic subspace. Furthermore there exists a natural 
grading operator on Tt, given by 

We have Q2 = iQi • X- Both Qi and Q2 anticommute with \ an d are therefore 
called supercharges of an abstract super symmetric quantum mechanics (cf. [lOj, 
sect. 8.3 or [TT] . [T2] respectively). We note that many of the susy properties 
follow already from the existence of this abstract structure. 

Remark 4.1 Note that x is both self adjoint and unitary with 

X = X* = X~ l , X 2 = 1 (35) 
The projectors on Tio and 7i\ are given by 

P = 1/2(1 + X ), Pi = 1/2(1 -x) (36) 
respectively, as for ordinary continuum Dirac operators. 
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5 Some Graph-Spectral Properties following from 
Susy 



We indicated already in [T] (see also [2]) that a variety of spectral properties on 
graphs are encoded in our Laplace or Dirac operator. For a finite, connected (for 
reasons of simplicity only) graph we have for example 

Observation 5.1 The following operator kernel, range properties hold. 

Dim (Rg d*) = n - 1 , Dim (Ker d*) = Dim(Rg rf) x ) = ^ vf - (n - 1) > 

(37) 

Dim(Rg d) = Dim(Ker d*) 1 ) = n - 1 (38) 

that is, Rgd* and Rgd have the same dimension, (n — 1), and the dimension 
of Kerd is one (n designates the order of the graph i.e. the number of vertices). 
Note that Y1 V T = J2 v i Ut = J2 v i/2 = ^(edges) is the dimension of Hi as each 
(directed) edge occurs as an ingoing edge for exactly one node. 

To give an idea how these results can be proved, take for example the first state- 
ment. Kerd is spanned by the vector x % as eacri edge, occurring as an in-edge 
for, say, Xi occurs as an out-edge for some other Xj, hence di^Xi) vanishes. 
With Rg A* = (Ker A) 1 - the result follows and the other results follow from sim- 
ple vector space mathematics and properties of the adjoint (note for example that 
for finite dimensions Dim(Ker d) + Dim(Rg d) = n)). 

For general infinite graphs our susy structure allows to infer more interest- 
ing spectral results. Note that, to keep matters simple, we restricted ourselves 
to graphs with globally bounded node degree. Therefore all our operators are 
bounded and there are hence no operator domain problems. In the following 
the polar decomposition of an operator turns out to be useful. With A a closed 
operator from a Hilbert space, Xi, to a Hilbert space, X2 we have the following 
representation 

A = S-\A\ (39) 

with \A\ = (A* A) 1 ' 2 a positive operator from Xi to Xi, S a partial isometry, 
mapping Rg | A | iso metrically onto Rg A (see for example JH] ) • 

We write d as d = S\d\ with \d\ = (d*d) 1/2 . We hence have d* = \d\S* and 
with dd* = \d*\ 2 = S\d\ 2 S* and consequently \d*\ = S\d\S* (the square root): 

d = S\d\ = \d*\S , d* = \d\S* = S*\d*\ (40) 

For our Q\ := D we have 

Lemma 5.2 The polar decomposition of Qi = D is 

fO d*\ ( S*\d*\\ fO S*\ f\d\ \ ,„ 1N 

Q ^={d o) = [s\d\ ) = {s 0H0 W ( ] 
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We furthermore have 
Lemma 5.3 

Ker Qi = Ker d © Ker d* = Ker Q 2 (42) 

with corresponding formulas holding for the respective orthogonal complements. 
For a finite connected graph we hence have (by observation ^. 1)) 

Dim{Ker Q lj2 ) = vf - (n - 1) + 1 = #(edges) - (n - 2) (43) 

For the proof note that for e.g. A bounded Ker(A*A) = Ker(A) as Ker(A*) = 
Rg(A)\ 

Proposition 5.4 On (Kerd) 1 - © (Kerd*) 1 - we have 

fd*d \ _ fS*dd*S \ 

^ dd*J ~ \ Sd*dS J 1 ' 

That is, d*d on Rg d* is unitarily equivalent to dd* on Rgd. This is a generalisa- 
tion to infinite dimensional spaces of a previous result concerning the dimension 
of the respective subspaces. It follows that the spectra of d*d and dd* coincide 
away from zero! We have 

d*df = Ef dd*Sf = ESf (45) 

Therefore the eigenvalue spectrum of Hs away from zero is at least twofold de- 
generate. 

H s (f,g) T = E{f,g) T => d*d f = E f and dd* g = E g (46) 

hence, (/, 0) T and (0, g) T are eigenvectors of Hs to the same eigenvalue and are 
at the same time eigenvectors of \ to the eigenvalues ±1. 

Furthermore, as Hs — Q\ — Q\ and {Qi,Q2} = 0, certain combinations of 
the above eigenvectors yield common eigenvectors of the pairs Hs, Qi- 

Proposition 5.5 The eigenvalues away from zero of Hs = Q\ are at least 
twofold degenerate. With E > being an eigenvalue of H s = Q\ = Q\ with 
(f,g) T the corresponding eigenvector of Q\, the eigenvalue of Qi is +E 1 ^ 2 or 
—E 1 ! 2 and, due to the anticommutation relation between Qi and Q2, Q2 {fi9) T 
is another eigenvector to Qi with eigenvalue —E 1 ! 2 or +E 1 ^ 2 . We thus see that 
both Qi and Q2 have a symmetric eigenvalue spectrum. 

We can make the result a little bit more explicit by making canonical choices. 
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Observation 5.6 With {f,g) T an eigenvector of Qi to eigenvalue A (and hence 
(f,0) T , {0,g) T "pure" eigenvectors of H$ to eigenvalue X 2 ), (f,—g) T is an eigen- 
vector to eigenvalue —A. Correspondingly, a straightforward calculation shows 
that (if, g) T is an eigenvector to Q2 with eigenvalue A and (if, —g) T the eigen- 
vector to eigenvalue —A. All these vectors are eigenvectors of Hs to eigenvalue 
E = A 2 . However, only the pairs belonging to Q\ or Q 2 , respectively, are linearly 
independent. We see that Q 12 necessarily mix the pure bosonic and fermionic 
eigenstates of H s (see also sect. 2.3.3 of TWp. 

Boilt down to the two operators d*d and dd* we have that with / being an 
eigenvector of d*d, df is an eigenvector of dd* to the same eigenvalue with a 
corresponding result for dd* , g and d*g (dd*df = dd* d f = Edf). Hence, all 
eigenvectors of dd* are of the form g — d f, f an eigenvector of d*d, both belonging 
to the same eigenvalue (see also ^T] sect. 5.2.3). Furthermore, the eigenvectors, 
(/, g) T of, for example, Q\ are characterized by the following symmetry property: 

d*g = \f , df = \g (47) 

We conclude from our preceding findings that the susy structure relates the spec- 
tral properties of d, d* , d*d, dd* to each other. On the other hand it seems to say 
(at least as far as we can see) not much about the spectrum of, for example, d*d 
as such. 

6 The Zero Eigenspace 

In supersymmetric quantum mechanics the zero eigenspace of Hs is particularly 
interesting and is associated with the notion of supersymmetry breaking (see for 
example ^D])- Our above lemma 15. HI makes an explicit statement about the 
dimension for a finite connected graph. We see that the eigenspace is in many 
cases highly degenerate and always has a dimension bigger than zero. This can 
be seen as follows. Each connected graph contains a spanning tree (see e.g. |14j). 
A finite tree has (n — 1) edges (see below). Hence each connected graph contains 
at least (n — 1) edges so that 

Dim(KerH s ) > 1 (48) 

In case of a finite tree we have the following. 

A finite tree of order n has (n — 1) edges. This can most easily be seen by 
choosing a base vertex, xq, and then starting from the outer vertices. Each edge 
corresponds to exactly one vertex, ending with the base vertex, which does not 
correspond to an edge. 

Lemma 6.1 For a finite tree the preceding formula reduces to 

Ker d* = (n - 1) - (n - 1) = (49) 
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We infer that Dim(Ker Hs) = 1. 

The situation is more complicated for infinite connected graphs. In this case 
the kernel of d is zero. The formula 

= df = J2(f*-fk)d lk (50) 

ik 

implies f) = for all pairs (i, k) occurring in the sum. The only normalisable 
vector in the infinite case has fa — for all %. 

On the other hand, for an oriented or (more generally) directed graph each 
cycle lies in the zero eigenspace of d*. More specifically, denoting the cycle by the 
edge sequence x^, . . . , x^ with x^ linked to x^ , we choose the following vector in 
Hi-. 

g ■= Y^[di„,iv+l] (51) 
v 

with [di^uz+i] denoting either d iv ^ v+ i if the edge is pointing from x iv to x iu+1 or 
—di v+1> iv if it points in the opposite direction. Applying d* to this vector yields the 
nullvector in Hq. Furthermore we state without proof that (genuinely) different 
cycles are linearly independent and span the kernel of d*. This can be shown by 
exploiting the existence of a spanning tree (see, for example [Hj, p. 53). We then 
have 

Observation 6.2 For an infinite connected directed graph the zero eigenspace of 
Qi, Q2 or Hs is purely fermionic and consists of the cycle space. 

This implies that the susy-Hamiltonian for an infinite tree has no non-trivial 
zero eigenvectors. We provide a separate proof for this statement as it employs a 
possibly useful technical property. We again pick a base vertex, Xq, and construct 
the spheres 

FiM := {xi , d{x , Xi) = 1} (52) 

with d denoting the canonical graph metric. We can infer that for a tree all the 
(x{ — Xk) so that linked by an edge are linearly independent in Ho and, 

together with xq span Ho- We then have that 

= d* g = ^ g ik {xj - x k ) (53) 

implies g^ = for all occurring pairs. 

Remark 6.3 This fact can also be exploited for general graphs by using a span- 
ning tree. 
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